CALCULATION OF TURBULENT SHEAR STRESS IN SUPERSONIC 


BOUNDARY LAYER FLOWS 
by 

Chen-Chih Sun 
and 

hL 

Morris E. Childs' 

University of Washington, Seattle, Washington 

Turbulent shear stress distributions for supersonic boundary layer 
flows have been computed, from experimental mean boundary layer data. The 
computations have been made by numerically integrating the time averaged 
continuity and streamwise momentum equations. The computational method 
is different from these previously reported in that integrated mass and 
momentum flux profiles and differentials of these integral quantities 
are used in the computations so that local evaluation of the streamwise 
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flows upstream and downstream of shock wave -boundary layer interactions 
and for both two-dimensional and axisymmetric flows. The computed 
results are compared with recently reported shear stress measurements 
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which were obtained by hot wire anemometer and laser ve lo c fmete r ‘ t e'qh 




niques . For most of the results shown reasonably good agr eeraerN^iis "found 
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between the computed and- the measured results. The computed values of 
shear stress are quite sensitive to rather small differences in the mean 
flow properties of the boundary layer, indicating that exceptionally 
accurate mean flow data are required if reliable shear stress distribu- 
tions are to be obtained from the mean data. Computations have also been 
made of eddy viscosity and mixing length distributions. These are 
compared with results based on the direct measurements of shear stress. 
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SYMBOLS 

a = a constant, see Eq. (13) 

C _ constant in Law of the Wall (usually equals 5.1) 

C.c = skin friction coefficient 

r 

K = a constant, see Eq. (13) 

k = a constant, see Eq . (18) 

£ = mixing length, see Eq. (17) 

M = Mach number 

F - pressure 

r - distance normal to centerline 

R = radius of the duct 

Ty = total stress tensor 

u = time averaged velocity in primary flow direction 

u* - (u /o 1/2 ) arc sin (o 1/2 u/u ) 

e g 

1 / 2 . 

u = friction velocity (t /p ) 1 

T w w 

v = time averaged velocity normal to centerline 

x = distance parallel to centerline • 

y = R-r 

Y ® ratio of specific heat (1,4 for present study) 

6 = boundary layer thickness 

e - eddy viscosity 

p = molecular viscosity 

n = y/6 

n = coefficient of wake function, see Eq. (13) 

p = time averaged density 



a = I ( y - 1 )/ 2]>£/{1 + Ky - 1 )/ 2 ] M *} 

t - shear stress 


Subscripts 

e - boundary layer edge condition 

w = wall condition 

» - free stream condition 

Superscript 

<( ) ! >- time averaged fluctuation value 
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X* INTRODUCTION 

In the study of supersonic turbulent boundary layer flow the 

turbulent shear stress distribution has always been of great importance 

and interest. The direct measurement of the turbulent shear stress is, 

however, quite difficult. A natural alternative is to compute the shear 

from experimental mean flow data by numerically integrating the momentum 

equation. Such computations have been performed in recent studies by 

1 2 3 

Bushnell and Morris , Horstman and Owen , and Sturek . Since directly 
measured data with which to compare the computed results were not avail- 
able, it was not possible to check the validity of the computations nor 
of simplifying assumptions made in connection with the computations. 

Recently, some turbulent shear stress data measured directly 

by using hot wire anemometer and laser velocimeter techniques have been 

4 5 

reported by Rose and Johnson * . The measurements were made upstream 
and downstream of an adiabatic unseparated interaction of an oblique 
shock wave with the turbulent boundary layer on the flat wall of a two- 
dimensional, M ro = 2.9 wind tunnel. The shock wave was generated by a 7° 
wedge. The turbulence data obtained from the two independent systems of 
measurement were in reasonably good agreement, indicating that the data 
should be reliable. 

In another study by Rose^ of the interaction of a conical shock 
wave with a turbulent boundary layer on the wall section of an 
axisymmetric wind tunnel, a hot wire anemometer has been used to make 
measurements of the turbulent shear stress and other turbulence quantities. 
The flow was adiabatic. For both of these investigations the turbulent 
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shear stress distributions at stations upstream and downstream of the 

shock wave-boundary layer interaction have been computed from the mean 

flow pitot pressure profiles obtained in the studies. Eddy viscosity 

and mixing length distributions have also been computed for the two 

flows. For the two-dimensional interaction investigated by Rose and 

45 

Johnson * , computed results show reasonably good agreement with the 
results obtained by direct measurement. Eor the axisymmetric flow 
studied by Rose^ computed and measured results show good agreement 
upstream of the interaction but differ considerably at the downstream 
station. 


II. BASIC EQUATIONS AND BOUNDARY CONDITIONS 

The time-averaged equations for the conservation of mass and 
momentum for steady compressible turbulent boundary layer flow in an 
axisymmetric channel are, respectively, 


3 / \ I 3 
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V XX 


( 3 ) 


T rx “ ^ T v ^rx ” <u,v,> + u < P , v , > + v <p T u’>) 

with t v representing the viscous stress. 


(A) 
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If we assume that [v <p*u f >{ << |p <u’v T> | and <<: I ;) x and 

transform to an x— y coordinate system the continuity and momentum equations 


may be combined and integrated in a direction normal to the surface to 
yield 
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where 
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( 6 ) 
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Rearrangement of equation (5 ) gives 
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Equation (8) becomes applicable to two-dimensional flow as R -* ®. 

The normal stress, T , which appears in equation (8), is not known 
from mean profile data. If its effect is to be considered in the computa- 
tions, it is necessary then to make assumptions regarding its magnitude. 

In earlier studies, and for most of the results to be discussed here, the 
streamwise gradient of has been neglected. However, by assigning 
arbitrary values to the derivative, its effect on the computed results 
may be examined. This has been done in one instance as will be discussed 
in the section on results. 

Knowledge is also required of the static pressure distribution in the 
boundary layer. In many studies of supersonic boundary layer flow no 
attempt is made to measure the static pressure variation normal to the 
wall, even though, for some adverse pressure gradient flows, the variation 
may, in fact, be rather large. In most instances the static pressure at 
the boundary layer edge may be determined with confidence. If this is done, 
the normal pressure variation may then be represented in approximate 
fashion by assuming a linear distribution between the wall static pressure 
and the pressure at the boundary layer edge, i.e., 
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P “ P w + n(P e - p w ) (9) 

Examination of equation (8) shows that an accurate value of boundary 
layer growth rate is very important for the calculation of the shear 
stress. However, precise determination of the boundary layer thickness 
from experimental mean data is difficult. It is even more difficult to 


evaluate the boundary layer growth rate accurately. This problem may be 
avoided by using the condition that the shear stress diminishes to zero 


at the boundary layer edge and solving equation (8) for d6/dx: 
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In equation (10) a linear pressure variation normal to the wall has been 
assumed. 

If the streamwise gradient of the normal stress is ignored and if the 
flow is assumed to be locally similar so that the nondimensional boundary 
layer properties u/u^ and p/p 0 are functions of n only, equation (8) 
reduces to 
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and equation (10) becomes 
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Equation (11) is identical to the expression used by Bushnell and Morris 
in their computations of turbulent shear stress. 

Before solving for the shear stress from either equation (8) or (11) 
it is also necessary to know the coefficient of skin friction. This may 
be obtained by using the wall-wake velocity profile proposed by Sun and 

n 

Childs . The wall-wake velocity profile for isoenergetic flow may be 
written as 
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As has been reported in Reference 7, the use of a - 1 in equations (13) and 
(14) results in a profile which provides a good representation of the 
boundary layer velocity distribution. The method of least-squares may be 
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used to fit the wall-wake profile to the experimental mean velocity pro- 
files to provides values of and to provide a smoothed representation 
of the mean velocity distribution. 

The eddy viscosity and mixing length distributions are also of great 
interest in the study of turbulent boundary layer flow. In terms of the 
shear stress as computed from mean flow data, the eddy viscosity may be 
written as 
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and the mixing length as 
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In the region close to the wall the mixing length may he approximated by 


Z = ky 


( 18 ) 


where k is a constant. 

In terms of turbulence measurements the eddy viscosity may be expressed as 
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and the mixing length as 
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A comparison of eddy viscosity and mixing length values as determined from 
turbulence measurements and from mean flow shear stress computations is 
given in the section on results. 


X 
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III. RESULTS 

A computer program has been developed to carry out the required 
numerical calculations. In carrying out the computations for the flow 
downstream of the shock wave-boundary layer interactions, the departure 
from local similarity has been taken into account. For purposes of 
comparison, however, computations based on local similarity have also been 
made. The effect of static pressure variation normal to the wall has also 
been considered for the downstream stations by assuming a linear variation 
in pressure. In these computations the static pressure at the boundary 
layer edge has been computed from the free stream total pressure and pitot 
pressure, with appropriate allowance made for the loss in total pressure 
across the shock system. As was mentioned in the previous section, the 
mean velocity profiles may be smoothed by using a least squares fit of the 
wall— wake velocity profile to the experimental profiles. Computations have 
been made for both smoothed and unsmoothed profiles. The integration 
process for the smoothed data should be more accurate because it is possible 
to use smaller step sizes. On the other hand, the shear stress computed 
from the smoothed profiles is valid only if the wall— wake profile provides 
a good representation of the actual velocity distribution. 

In obtaining the density and velocity profiles from the pitot pressure 

/ C 

profiles for the investigation by Rose and Johnson * , the total tempera- 
ture was assumed to be constant across the boundary layer. In the axisym- 

g 

metric study by Rose total temperature measurements were available. For 
this case calculations of the shear stress distributions were made for both 
the measured total temperature distribution and for a constant total tem- 
perature. The results differed from each other only slightly. 
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Figure 1 shows shear stress distributions computed for an upstream,. 

(x - 5.375 cm) station in Johnson's and Rose’s investigation, along with 
their shear stress data from the hot-wire anemometer and the laser 
velocimeter measurements. The computed shear values which are shown 
have been obtained from equation (11) and auxiliary equation (12), i.e,, 
under the assumption of local similarity. As is shown the calculated 
results agree quite well over much of the boundary layer, with the meas- 
ured results obtained with the laser velocimeter. The differences between 
the calculated results and the hot-wire results are greater. For both 
the hot-wire and laser velocimeter measurements the peak values of shear 
stress are seen to occur substantially farther from the wall than is 
observed for the calculated distributions. Also shown in the figure are 
values of the wall shear stress as determined by a least squares fit of 
the wall-wake profile to the mean data and as measured by a Preston tube. 
The agreement between the two shear stress values is good. 

The wall-wake profile as given by equation (13) can be reduced to 

, g 

a profile proposed earlier by MaiseO and McDonald if II is taken to be 
equal to 5.0 and if a «>. For purposes of: comparison, the Maise- 
McDonald profile has been used to compute shear stress values. As is 
shown, this results in shear values which are substantially lower than 
those obtained with a = 1 in equation (13) or with the unsmoothed data. 

Computations have also been carried out for the data of Rose and 
Johnson for a station downstream of the shock wave-boundary layer inter- 
action (x = 9,375 cm). Several sets of computations have been made from 
the mean data. In one set local similarity has been assumed. Equations 
(11) and (12) have then been used to carry out computations for both 
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wall-wake smoothed profile and a profile as measured. In. other computations 
local similarity of profiles has not been assumed. For these cases re- 
sults have been obtained for both smoothed and unsmoothed profiles and both 
linear and constant static pressure distribution across the boundary layer. 
The results are shown in Figures 2 and 3. Figure 2 shows that the meas- 
ured shear stress values obtained by the two experimental techniques are 
in quite good agreement. These in turn agree reasonably well with 
the computed values obtained by using smoothed nonsimilar profiles and the 
assumption of a linear pressure variation across the boundary layer. There 
is considerable difference between the results for the smoothed and un- 
smoothed profiles with the peak shear stress value for the unsmoothed pro- 
file occurring lower in the boundary layer. As was stated earlier, the 
numerical integration process should be more accurate for the smoothed 
profiles. On the other hand, the results obtained from the wall-wake 
profiles are valid only if the profiles provide accurate representations 
of the actual velocity distributions. A comparison of the computed 
velocity distributions with the corresponding wall-wake profiles is shown 
in Figure 4. The differences between the two are small, pointing to the 
sensitivity of the computed shear stress values to small velocity differ- 
ences. 

As is shown in Figure 3, the results obtained under the assumption 
of local similarity are markedly different from those determined when 
similarity is not assumed. The peak shear stress levels computed assuming 
local similarity are less than half the values computed when similarity 
is not assumed. Furthermore, the shapes of the shear stress distribution 
curves are quite different. These results occur even though the differences 
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between profiles at two closely spaced stations are small. Figure 5 
shows comparisons of density and velocity profiles from Rose’s and 
Johnson’s study for streamwise stations located approximately one boundary 
layer thickness apart. Although the profiles appear at first glance to be 
quite similar, striking differences are found for the computed shear 
stress distributions. 

g 

In the investigation by Rose of an adiabatic axisymmetric shock 
wave-boundary layer interaction the free stream flow was nominally at M - 4 
and the shock wave was generated by a 9-degree half-angle cone placed at 
zero angle of attack on the centerline of the 5.28 cm diameter tunnel. The 
shock strength was just below that required to produce boundary layer 
separation. For this flow a substantial adverse pressure gradient existed 
downstream of the interaction. As was stated earlier, the turbulence data 
were obtained with a hot-wire anemometer. 

Figure 6 shows results for a station upstream of the interaction 

where there was no pressure gradient. The shear stress distribution has 

been computed from mean velocity and density profiles under the assumption 

of local similarity and constant static pressure across the boundary layer. 

As is shown, reasonably good agreement between computed and measured shear 

stress values is obtained, although the differences observed for the Rose- 
A 5 

Johnson * data are found to exist here as well. That is, the hot-wire 
anemometer results show higher levels than the computed values in the outer 
part of the boundary layer, and the peak value of the computed shear stress 
is much closer to the wall than was found to be the case from the direct 


measurements. 
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Figure 7 shows results for a station (at x = 10.16 cm) downstream of 
the interaction. As was the case with the Rose- Johnson investigation, 
several sets of computations have been made from the mean data. For all 
the results shown the assumption of a linear static pressure variation 
across the boundary layer has been assumed. Results have been obtained 
for both smoothed and un smoothed profiles, and computations have been 
carried out for both local similarity and nonsimilarity of velocity and 
density profiles. In addition, allowance has been made in one set of 
calculations for the strearawise gradient of the normal stress, 9T /3x. 
Johnson and Rose found in their study of the two-dimensional interaction 
that this gradient was small and no allowance was made for it in the compu- 
tations. However, for the axisymmetric flow Rose reported that stream- 
wise normal stress gradient reached levels as high as 30 percent of the 
local streamwise pressure gradient. The distribution of normal stress 
gradient as reported by Rose may be represented approximately by 
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n , 3P .2 
0.3 sxn* r}7i 


( 21 ) 


Computations have been made assuming the distribution given by equation 

( 21 ) . 

As is shown in Figure 7, the difference between the measured and 
computed shear stress values is substantial. The peak value of the meas- 
ured shear stress is about twice the peak value computed when local simi- 
larity is not assumed. As was true for the downstream station in the Rose- 
Johnson study, peak levels of the computed shear stress are seen to occur 
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considerably closer to the wall than is observed in the measurements. It 
has not yet been possible to determine the causes for the rather large 
differences between the peak levels of the computed and measured shear 
stress distributions. Additional hot-wire studies are underway for the 
axisymmetric flow, as are attempts to improve upon the computational 
techniques. 

4,5 

As was the case for the downstream station in the Rose-Johnson 
investigation, the assumption of local similarity leads to shear stress 
levels which are on the order of only one-half those determined when 
similarity is not assumed. As is shown in figure 7, consideration of the 
streamwise normal stress gradient in the computations results in shear 
stress levels which differ only very slightly from those obtained when 
the gradient is ignored. 

£ 

It should be pointed out that Rose , in obtaining his downstream- 
data, moved the cone along the axis of the tunnel and obtained the pitot 
measurements at a fixed axial station. The question arises, then, of 
how the computed shear stress distributions would compare with those 
obtained if the cone had been kept at one station and the pitot tube had 
been traversed axially. As a check on this, data were examined from a 
study by Seebaugh^ in which the cone was held fixed and the probe was 
translated. Seebaugh’s study was conducted in a 5.16 cm diameter round 
tunnel with a free-stream Mach number of 3.78. He used a 10-degree half- 
angle cone. The results are shown in Figure 8. As is apparent, the shear 
stress distributions are quite comparable to those obtained from Rose’s 


mean flow data. 
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The distributions of eddy viscosity and mixing length for the inves- 
tigations under consideration are shown in Figures 9-12. The results, 
of course, reflect the differences between the measured and computed 
values of shear stress as discussed earlier. Computed results are shown 
for both smoothed and unsmoothed profiles. For the stations upstream of 
the shock wave-boundary la} 7 er interactions the shear stress distributions 
used in obtaining the values of eddy viscosity and mixing length have 
been determined under the assumption of local similarity. For the down- 
stream stations, the lack. of local similarity in the velocity and 
density profiles has been taken into account and a linear static pressure 
variation across the boundary layer has been assumed. 

The eddy viscosity distributions for the Rose-Johnson investigation 
are shown in Figure 9. Near the wall moderately good agreement is 
observed between eddy viscosities based on computed shear stress values 
and those based on measurements. From y/6 ~ 0.2 outward., however, the 
differences become substantial. The values based on experimental shear 
stress show considerable scatter as do those based on the point-by-point 
mean data. In contrast the results based on the wall-wake smoothed 
velocity profile are smooth, as would be expected. 

4 5 

The mixing length distributions for the Rose-Johnson ’ investigation 
are shown - in Figure 10. Near the wall the values of mixing length are 
described reasonably well by the relationship Z ~ 0.4y. From y/6 - 0.2 
outward a certain amount of scatter in the results is apparent. For the 
upstream station (x - 5.375 cm), however, the computed distribution based 
on the wall-wake profile agrees with the results based on turbulence meas- 
urements reasonably well. For the downstream station (x ~ 9.375 cm) the 
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values of f/6 based on the wall-wake profile are consistently higher than 
those based on the turbulence measurements or on the unsmoothed mean dat«.. 

The eddy viscosity distributions for the axisymmetrie flows studied 
by Rose 6 and by Seebaugh ^ are shown in Figure 11. As is shown the 
results based on the shear stresses computed from the mean flow data of 
Rose and Seebaugh are in good agreement. The values based on Rose’s 
hot-wire anemometer measurements are much higher, consistent with the 
higher shear stresses which he reported. 

Mixing length distributions for the downstream stations in Rose's 

and Seebaugh’ s investigations are shown in Figure 12, along with the 

results based on Rose's hot-wire anemometer measurements. Again the 

results agree well with the expression & — 0.4y for the region near the 

3 

wall. This is in contrast to the result reported by Sturek who concluded 
for adverse pressure gradient flow £ = 0.65y appeared to fit his data near 
the wall. Sturek' s studies were conducted for flow, along a compression 
surface. The computed values of mixing length, £/6, in the plateau 
region are considerably different for the two-dimensional flow and the 
axisymmetrie flows. Upstream of the shock wave interaction, the value is 
about 0.07 for axisymmetrie flow while the value is about 0.12 for the 
two-dimensional flow. The computed values downstream of the interaction 
are somewhat higher than those upstream for both the two-dimensional and 


axisymmetrie flows. 
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IV. CONCLUSIONS 

A method of computing shear stress distribution from experimental 
mean profile data in compressible turbulent boundary layer flow has been 
developed. The method is different from those previously reported in 
that integrated mass and momentum flux profiles and differentials of 
these integral quantities are used in the computations so that local 
evaluation of the streamwise velocity gradient is not necessary. The 
method has been found to yield results which are in reasonably good agree- 
ment with directly measured turbulence data for two-dimensional adiabatic 
boundary layer flow in the regions upstream and downstream of an oblique 
shock wave interaction. The computed results are very sensitive to the 
accuracy of the numerical integrations required in the computational 
procedure and to the mean property distributions in the boundary layer. 

The assumption, of local similarity may cause large errors in computed 
shear stress values for flows subjected to pressure gradients, even though 
adjacent profiles of the mean properties appear to be quite similar. 

The shear stress levels are quite sensitive to the static pressure dis- 
tribution normal to the wall. The effect of the streamwise gradient of the 
normal stress on the computed results is small and apparently may be 
neglected. The value of the constant k in the expression £ = ky for 
the mixing length in the region near the wall remains close to 0.4 for 
adverse pressure gradient flows along a flat surface. In view of the 
rather substantial differences between results based on mean flow measure- 
ments and those based on turbulence measurements in an axisymmetric 
adverse pressure gradient flow, further study of this flow is needed. 
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Figure 1. Turbulent shear stress distribution upstream 
of a shock wave-boundary li^yer interaction, 
two-dimensional tunnel [4 ] 


Figure 2. Turbulent shear stress distribution 
downstream of a shock wave-boundary 
layer interaction, two-dimensional 
tunnel [5] 


Figure 3. Effect of static pressure distribution and 

self-similarity on shear stress distribution 
in a turbulent boundary layer downstream of 
a shock wave interaction [5] 


Figure 4. Velocity profiles and wall-wake representations 
downstream of a shock wave-boundary layer 
interaction [5] 


Figure 5. 


Density and velocity profiles downstream of a 
shock wave — boundary layer interaction rsi 


Figure 6*. Turbulent shear stress distribution upstream 
of a shock wave-boundary layer interaction, 
axisymmetric tunnel [6] 


Figure 7. Turbulent shear stress distribution downstream 
of a shock wave- boundary layer interaction, 
axis ymme trie tunnel [63 


Figure 8. Turbulent shear stress distribution downstream 
of a shock wave-boundary layer interaction, 
axisymmetric tunnel [9] 


Eddy viscosity distribution for the two-dimensional 
flow of Rose and Johnson [4,5] 


Figure 9. 
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Figure 10. Mixing length distribution for the two-dimensional 
flow of Rose and Johnson [4,5] 


Figure 11. Eddy viscosity distributions for the 
axisymmetric flows ox Rose [6} and 
Seebaugh [9] 


Figure 12, 


Mixing length distributions for the 
axis ymme trie flows of Rose [6] and 
Seebaugh [9] 
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Measured - [- p<u'v ? >/p u^] x 10 3 
O laser velocimeter [4] 

A - hot-wire anemometer [4] 
Computed - [(x - p3u/3y)/p u^] x 10 3 
D mean data, eq. (11) 


wall-wake presentation of mean 
data, eq. (11) 


Mai s e-McDonald profile repre- 

sentation of mean data [8] 

q (^}- T w computed from wall-wake profile 
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Computed - f (t - y3u/3y)/p e u^] x 10 3 

Q mean data, linear static pressure distribution 
Eq. (8) 

wall-wake representation of mean data, linear 

static pressure distribution, Eq. '(8) 
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A hot-wire anemometer [6] 
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O mean data, Eq. (11) 

wall-wake representation of mean data, Eq. (11) 

(T) t computed from wall-wake profile 
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wall-wake representation of mean data, Eq. (8) 



0.0 



cm 


3 



( 8 ) 

( 11 ) 


* i 



(e/u 6) x 10 3 


33 


Eddy Viscosity based on measured shear stress 
# x = 5.375 cm, laser velocimeter [4] 

O x * 9.375 cm, laser velocimeter [5] 

A X s * 5.375 cm, hot-wire anemometer 

A x = 9.375 cm, hot-wire anemometer 

Eddy Viscosity based on computed shear stress 

G x = 5.375 cm, computed from mean data, Eq. (8) 

□ x - 9.375 cm, computed from mean data, Eq. (8) 

x - 5.37S cm, computed from wall-wake representation of 

mean data, Eq. (8) 

x = 9.375 cm, computed from wall-wake representation 



- i 



&¥■ 


Mixing Length based on measured shear stress 
O x = 5.375 cm, laser velocimeter [4] 

O x - 9.375 cm, laser velocimeter [5] 

A x * 5.375 cm, hot-wire anemometer 

A x - 9.375 cm, hot-wire anemometer 

Mixing Length based on computed shear stress 

D x = 5.375 cm, computed from mean data, Eq. (8) 

Q x = 9.375 cm, computed from mean data, Eq. (8) 


x = 5.375 cm, computed from wall-wake representation 
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Eddy Viscosity based on measured shear stress 
A x ■ 6.60 cm, hot-wire anemometer [6] 

A x = 10.16 cm, hot-wire anemometer [6] 

Eddy Viscosity based on computed shear stress 

0 x = 6.60 cm, computed from mean data [6] 

x = 6.60 cm, computed from wall-wake representation 

of mean data 

O x = 10.16 cm, computed from mean data [6] 

x = 10.16 cm, computed from wall-wake representation 

of mean data [6] 

CD x = 9.62 cm, computed from mean data [9] 
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Mixing Length based on measured shear stress 
A x = 6.60 cm, hot-wire anemometer [6] 

A x “ 10.60 cm, hot-wire anemometer [6] 

Mixing Length based on computed shear stress 

# x = 6.60 cm, computed from mean data [6] 

x = 6.60 cm, computed from wall-wake representation 

of mean data [6] 

O x “ 10.16 cm, computed from mean data [6] 

x = 10.16 cm, computed from wall-wake representation 

of mean data [63 

□ x - 9.62 cm, computed from mean data [9] 
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